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The 7a discrepancy between the proton rms charge radius from muonic hydrogen and the 
CODATA-f value from hydrogen spectroscopy and electron-scattering has caused consider- 
able discussions. Here, we review the theory of the 2S-2P Lamb shift and 2S hyperfine split- 
ting in muonic hydrogen combining the published contributions and theoretical approaches. 
The prediction of these quantities is necessary for the determination of both proton charge 
and Zemach radii from the two 2S-2P transition frequencies measured in muonic hydrogen 



I. INTRODUCTION 

The study of the energy levels in hydrogen-like atoms (hydrogen, muonium, positronium) as 
well as free and bound g-factors provides the most precise verification of quantum electrodynamics 
(QED). 

Comparison of the measured transition frequencies in hydrogen with theory is limited by the 
uncertainty of the proton structure, mainly by the uncertainty of its root mean square (rms) charge 
radius defined as r| = J d 3 rr 2 p(r) where p is the normalized charge density of the proton. Testing 
of QED in muonium and positronium does not suffer from uncertainties of the nuclear structure, 
but the ultimate experimental precision is limited by the finite lifetime of these systems. A way out 
to improve the test of hydrogen energy levels is given by spectroscopy of muonic hydrogen (/up), 
an atom formed by a muon and a proton. This provides a precise determination of te- 

Recently the comparison of the measured 2P^ 1 — 2S^° transition in muonic hydrogen j3| 

A£ 2 c ^ 2 _ 2gf ^ = 206.2949(32) meV (1) 

with the theoretical prediction based on bound-state QED 0-0] 

AE& F=a QQF=1 = 209.9779(49) - 5.2262 r 2 + 0.03474 meV (2) 

Zr 3/2 ~ /S l/2 V ' 

yielded te = 0.84184(67) fm Q]. Yet this value, although being an order of magnitude more precise, 
is 4% (or 7 standard deviations) smaller than the CODATA-10 value r g ODATA = 0.8775(51) fm 0] 
based on hydrogen spectroscopy and electron-proton scattering. 

A large effort has been devoted to solve the proton size puzzle ranging from studies of physics 



beyond the standard model [7H16I] to new measurements and reanalysis of electron-proton scattering 
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data [17H22j , from theories describing the proton structure at low energies |23l-l27l] , to onsets of new 
determinations of the Rydberg constant (Roo) at MPQ, LKB, NPL and NIST. Several studies have 
been dealing also with the theoretical prediction of the n = 2 energy levels in muonic hydrogen. 
The aim of this paper is to summarize the various contributions of these investigations updating the 
theoretical prediction of the muonic hydrogen 2S-2P Lamb shift and 2S-HFS. It can be anticipated 
here that no big error or additional contributions have been found which could solve the observed 
discrepancy between experiment and theory of the 2P-2S energy difference: 

discrepancy = AE™ P F =2 F=1 - AE^ 2 _ 2SF=1 (rg ODATA ) = 0.31 meV (3) 

Zr 3/2 ~ Z0 l/2 Zi 3/2 zo l/2 

corresponding to a relative discrepancy of 0.15 %. 

Experimentally two transition frequencies in muonic hydrogen have been measured, one starting 
from the 2S-triplet state v t = 2P|^ 2 — 2S^ 1 Q, 0], and the other from the 2S-singlet state 
u s = 2P|^ 1 — 2Sf^° [l| as shown in Fig. [TJ As presented in Sec. [TT] we can deduce from these 
measurements both the "pure" Lamb shift (AE^ = A£ , 2P 1/2 -2S 1/2 ) and the 2S-HFS splitting 
(AE'hfs) independently of each other. Comparing the experimentally determined AE-^ with its 
theoretical prediction given in Sec. IIIII yields an improved te value free of uncertainty of the 2S- 
HFS. Similarly comparing the experimentally determined A£hfs with its theoretical prediction 
given in Sec. IIVI results in the determination of the Zemach radius rz- 

Perturbation theory is used to calculate the various corrections to the energy levels involving an 
expansion of both operators and wave functions. The radiative (QED) corrections are given by an 
expansion in a, binding effects and relativistic effects in (Za), and recoil corrections in the ratio 
of the masses of the two-body system (m/M). Z = 1 is the atomic charge number and a the fine 
structure constant. The contributions related with the proton structure are in part described by 
an expansion in powers of te and rz- The book-keeping of these several corrections contributing 
to the /j,p Lamb shift and 2S-HFS is challenging because: 

• All corrections are mixed-up as a x (Za) y (m/M) z . 

• There are large finite-size and recoil (m/M ~ 1/9) corrections. 

• One cannot develop the calculation in a systematic way like in g — 2 for free particles. 

• Widely different scales are involved: the masses, the three-momenta and the kinetic energies 
of the constituents. 

• Different authors use different terminologies for identical terms. 

• Different methods are being used: Schrodinger term + Breit corrections versus Dirac term, 
Grotch- versus Breit-type recoil corrections, all-order versus perturbative, non-relativistic 
QED (NRQED) etc. 

In this study we summarize all known terms included in the Lamb shift and the 2S-HFS pre- 
dictions which are used in [l| to determine the proton charge radius and the Zemach radius. The 
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FIG. 1: 2S and 2P energy levels. The measured transitions i>t 01 and v 8 [l| are indicated together with 
Lamb shift, fine and hypcrfine splittings and finite size effects. The main figure is drawn to scale. The insets 
zooms in on the 2P states. Here, the mixing of the 2P(F=1) levels shifts them by ±<5 (see Eq.©). 



majority of these terms can be found in the works of Pachucki 0,0] , 
Karshenboim et al. 



Borie 



311 ] and Martynenko [32], 133J] . In addition Indelicato [3jj checked and improved 



Jentschura 



29 



301, 



many of the relevant terms by performing numerical integration of the Dirac equation with finite- 
size Coulomb and Uehling potentials. A detailed description of existing terms can be found in the 
review of Eides et al. [35 1. 



II. THE EXPERIMENTAL LAMB SHIFT AND 2S-HFS 



The 2S and 2P energy levels in muonic hydrogen are presented in Fig. [TJ They are given by the 
Bohr structure corrected for relativistic, hyperfme, radiative, recoil, and nuclear structure effects. 
The Lamb shift is dominated by the one-loop electron-positron vacuum polarization of 205 meV. 
This is much smaller than the n = 2 binding energy of 632 eV. Because we are measuring the 
Lamb shift between n = 2 states, the main term of the binding energy (Bohr structure) drops out 
and thus our results do not depend on the Rydberg constant. The experimental uncertainty of the 
Lamb shift is of relative order u r ~ 10~ 5 . Thus the various contributions to the Lamb shift have 
to be calculated to better than ~ 0.001 meV. 
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The two measured transition frequencies shown in Fig. [T] are given by 



hu t =£(2Pf 7 = 2 ) - E(2S{ 



1/2 > 



~-AE L + AE FS + -AE- 



2P3/2 

HFS 



1 



■A^HFS, 



(4) 



^ s =E(2Pf / = 1 )-S(2Sf / =°) 



,2P, 



--AE L + AE FS - ^AE" HF f + 6 + -AE HFS , 

with h being the Planck constant and AEfs = Ai?2P 3/2 ~~ A£2P 1/2 the fine splitting of the 2P-state. 
For the 2P fine structure we use the value 



33] 



A£ FS = 8.352082 meV, 



(5) 



to be compared with the values 8.3521 meV [28| and 8.351988 - 0.000052 r| = 8.351944 meV from 
Eq. (119) in [341 ]. The 2P 3 / 2 hyperfine structure is [33] 



A p 2F>3 / 2 

HFS 



3.392588 meV, 



(6) 



in agreement with 3.392511 meV (Indelicato and Mohr, 2010, unpublished). The 2P|^ 2 1 level is 
shifted upwards due to state mixing |3g] of the 2P(F=1) levels by [33] 



0.14456 meV. 



(7) 



The values given in Table 8 of [281 ] include this shift and deviate by less than 0.0002 meV from the 
values in Eqs. © and (0). 

A£l and A£hfs can be deduced independently of each other by forming the following linear 
combinations: 



-hv s + ^hu t =AE L + A^fs + ^A^fs 2 + \ 5 



3. . „ . „ 1 AF 2P ^ /9 1 

8 AjBh: 

~-AE L + 8.8123(2) meV, 



(8) 



hv s - hu t =AS H FS - A-Ehfs 2 + 5 



--AE UFS - 3.2480(2) meV, 



where the uncertainties of the constant terms correspond to uncalculated higher-order QED terms, 
and differences between the various authors. These uncertainties are negligibly small compared to 
the uncertainties of the measured transition frequencies [it and the uncertainties of the theoretical 
prediction for AE^ and AE^hfs ( see below). Finite-size effects for the 2P state are small but have 
been included in the theoretical prediction of both the fine and hyperfine contributions. For the 
fine splitting they have been computed perturbatively to be — 0.0000519 r 2 . = —0.00004 meV 28] 



and with an all-order approach -0.0000521 r| = -0.00004 meV |3J]. The finit e-size contributions 
to the 2P3/2 hyperfine splitting are < 10 -5 meV 34j. Hence the uncertainties of the constant terms 



in Eqs. (|8|) related with the proton finite-size are totally negligible. 
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III. THE LAMB SHIFT PREDICTION AND THE PROTON CHARGE RADIUS 



The main contribution to the Lamb shift in muonic hydrogen contrarily to hydrogen is given 
by the one-loop electron-positron vacuum polarization (eVP). The second largest term is related 
with the finite extension of the proton charge and is proportional to m^r^, where m r ~ 186m e is 
the reduced mass of the system. The third largest contribution is given by the Kallen-Sabry (two- 
loop eVP) diagrams and the forth by the muon one- loop self-energy summed with the one-loop 
muon-antimuon vacuum polarization (//VP). In has to be stressed that the observed discrepancy is 
larger than any other contribution besides the four explicitly listed here. Therefore a solution of the 
proton radius puzzle in the context of muonic hydrogen theory could arise either from one-loop VP, 
or proton structure effects, or fundamental problems in bound-state QED, but not from a missing 
or wrong higher-order contribution. Nevertheless it is necessary to compute all these higher order 
contributions to better than ~ 0.001 meV in order to exploit the accuracy of the measurements 
and to determine te with a relative accuracy of u r ~ 5 X 10 -4 corresponding to 0.0004 fm. The 
radiative, relativistic, binding and recoil corrections to AE^ have been summarized in Table [J 
whereas the proton-structure dependent contribution are summarized in Table [III 



A. Proton structure independent contributions 



The corrections to the energy levels predicted by the Schrodinger equation solved for the 
Coulomb potential are usually calculated using perturbation theory. Relativistic corrections are 
obtained from the two-body Breit-Hamiltonian. To check the validity of the perturbative approach 
the largest contribution (one-loop eVP) has been recalculated by integrating numerically the Dirac 
equation [2a . |34j, 1371 . |38j. The inclusion of the Uehling to the Coulomb potential yields the rela- 



tivistic all-order one- loop eVP correction given by the sum of #3 and #5 in Table HI This compares 
well with the perturbative results which correct #1 for relativistic and eVP iterations given by #2 
and #5. The entry #5 is the energy shift in second order perturbation theory related with the 
wave function distortion caused by one-loop eVP. The entry #7 represents the effect of the wave 
function distortion caused by one- loop eVP on the Kallen-Sabry contributions. 

The relativistic eVP corrections with account of the reduced mass dependence of 0.018759 meV 



was given m [30J. In [3l|] it was shown why it differs from the previous result of 0.0169 meV 
presented in #2. In #19 we give the difference between the full expansion in (m/M) of the 

(rcl) /2p 



1/2 



relativistic eVP correction (Eq. (4) in 3l|) and its lowest order (Eq. (6) in 3l|]): Ey P 
2S 1/2 ) - £^(2P 1/2 - 2S 1/2 ) = 0.018759 - 0.020843 = -0.002084 meV. Using these corrected 
values, the various approaches for the one- loop eVP contributions with relativistic-recoil corrections 
are consistent: starting from the Schrodinger equation Pachucki becomes 205.0074 + 0.018759 = 
205.0262 meV and Borie starting from the Dirac equation 205.0282 - 0.002084 = 205.0261 meV. 

For #11, the muon self-e nerg y correction to eVP of order a 2 (Za) 4 , there was an improvement 
summarized by Eq. (27) of [39J which includes contributions beyond the logarithmic term with 
modification of the Bethe logarithm to the Uehling potential. 

The term #12, eVP loop in self-energy, is part of #21 as can be seen from Fig. 22 in 40J]. Thus 



it was doubly-counted in "Old-Nature" compilation. 



B. Proton structure dependent contributions 

We first present the background required to understand the individual proton structure depen- 
dent contributions summarized in Table HU The main contribution due to finite nuclear size has 
been given analytically to order (Za) 6 by Friar 50|]. The main result is 28] 



2ira 



finite size 



*(o)| 5 



(r 2 ) --^m r Rf 2) + {Zaf (F REL + m 2 F N REL) 



(9) 



with 
-Frel 
FnreL 



(r 2 ) 



35 

7 - — + \n(Za) + (ln(m r r)} 
lb 



l (r 3 )( I )+/ REL + / REL 
o V 
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G 



+ ln(Za) 



+ — (r 2 )(r 2 ln(m r r)} 



tl + (r 3 )(r) + 1 (r 5 )( I } + jNREL + jNREL 

(10) 



where (r n ) = j d 3 r r" p(r) is the n th moment of the charge distribution, r E = \/ (r 2 ) ~ 0.84 fm 
thus the rms charge radius, R%\ the third Zemach moment (see below), m r the reduced mass, 7 
the Euler's constant and I^ Eh , jNREL^ /REL anc j jREL integrals which depend on the charge 
distributions. The terms proportional to (r 2 ) can be summarized for the 2S-state in muonic 
hydrogen as (b a + + b c + bd + b e )(r 2 } following the notation given in Table B.14 of 28f] . The 
terms proportional to bi„ bd and b e account for radiative finite-size corrections. 

The coefficient of the first term in Eq. ([9]) which is the leading finite-size effect, is for the 
2S state b a = 5.1973 meV/fm 2 [28(] calculated using the Schrodinger wave function. Relativistic 
corrections to finite-size are accounted for in -Frel- The coefficient proportional to (r 2 ) in .Frel is 
b c = 0.00181 meV/fm 2 . 

The second largest finite-size effect is related with the third Zemach moment 



Rf 2) = J d 3 r J dV r 3 PE (r-r r )p E {r') 



(11) 



which can be determined experimentally from the electric Sachs form factor Ge measured in elastic 
electron-proton scattering as 51] 

1 



R h 



f / d f,[GW)-l + - q \r 2 )]. (12) 
The resulting values are i?L = 2.71(13) fm 3 0] and J?L = 2.85(8) fm 3 Q. The second ter 



Eq. 



amounts to 



AF tnird Zemach — 0.0091i? 3 



(2)- 



m m 



(13) 



It is customary to express the third Zemach moment using the second moment of the charge 
distribution as 



f(r Z ) 



2\3/2 



(14) 
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where / is a constant which depends on the model for the shape of the proton. For an exponential 
charge distribution (dipole form factor) / = 3.79 is obtained. This has to be compared with 
/ = 3.47 for a Gaussian distribution to illustrate the sensitivity on the proton shape, and with the 
values extracted from the measured electric form factors / = 3.78(31) 51[ and / = 4.18(13) 18|. 
From Eq. ([13]) and Eq. {TI]) we get 



AEthird Zemach = 0.0365(18) (r 2 f' 2 



(15) 



(see first column in Table [TT]) where / = 4.0(2) [28j was used to account for the spread of the 
various values measured in scattering experiments. 

A solution of the proton radius puzzle assuming a large tail of the proton charge distribution 
result ing in an extremely large value 



data 



18 



20 



52 



531 ] has been ruled out by electron-proton scattering 



511 ] and by chiral perturbation theory [27j, |54j . 



The third Zemach contribution may be seen simplistically as a second-order correction in per- 
turbation theory given by the modification of the wave function caused by the finite size. In a 

0,0, 0,0]. The two- 



quantum field framework it is part of the two-photon exchange diagrams 
photon contributions to the finite-size are usually divided into an elastic part given by the R%\ 
term in Eq. ([9]), where the intermediate virtual proton remains on-shell, and an inelastic part (pro- 
ton polarizability contribution AE po \) where the virtual proton is off-shell. A unified treatment 
of such contributions can be achieved in modern quantum field theory using the doubly-virtual 
Compton amplitude which can be related to measured form factors and spin averaged structure 
functions using dispersion relations. Part of a subtraction term needed to remove a divergence 
in one Compton amplitude is usually approximated using the one-photon on-shell form factor 



A possible large uncertainty related with this approximation [25j has been recently discarde d by 



model independent calculations in the framework of heavy baryon chiral perturbation theory j^f . 

The total two-photon exchange contribution with finite-size amounts to A£tpe 

0.0333(20) meV fl, 



29 



551 ] . It results from the sum of an elastic part AE, 



cl 



0.0295(16) meV |23j, [55j, a non-pole term AE np = —0.0048 meV [23j, |55j, a subtraction term 
Ai? su b = —0.0041(10) meV [23J and an inelastic (proton polarizability) contribution AE po \ = 
0.0127(5) meV 55J]. Additional contributions from new models for the "off-shell" form factor have 
been limited to < 0.001 meV by quasi-elastic electron scattering 561 ] . The uncertainty of AE C \ 



accounts for various measured form factors (Kelly, AMT and Mainz 2010) 
The sum 

AE ei + AE np + AE snh = 0.0206(19) meV 



23 



(16) 



should be compared with AE'tMrd Zemach corrected for recoil corrections [29J which reduce the third 
Zemach contribution by a factor (0.018 meV)/(0.021 meV) 0] to 



AE, 



third Zcmach+rccoil 



0.018 
0.021 



• 0.0091^ 2) = 0.0217(15) meV 



(17) 



with the uncertainty arising from recoil corrections (5%) and the spreac. of the R%\ values 
The approximate scaling of 0.018/0.021 has been discussed by Jentschura 



291 ] based on Pachucki 
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to account for recoil corrections which are automatically included in the unified treatment of 
the two-photon exchange contribution AE'tpe- This recoil correction to two-photon finite-size 
contributions is ~ —0.003 meV, much less and of opposite sign compared to the term of 0.013 meV 
(#18 in Table Hi]) given by Borie following Friar. Comparison between Eq. (|16p and Eq. ()17|) shows 
that there is fair agreement between the results from the analytical formula of Eq. ([9]) and the 
quantum field approach via Compton scattering and dispersion relations. 

The sum of all the other terms in Eq. ([9]), which have been calculated assuming an exponential 



charge distribution and rE=0.875 fm, affects the Lamb shift by AE fi °j t ^°^ = —0.000123 meV [28fl . 



Even though the assumption of an exponential distribution may not be completely realistic for 
these higher-order contributions, this is sufficient [28j]. The smallness of this term may be qualita- 
tively understood in a perturbative-framework from Eq. ([9]) where higher moments of the charge 
distribution appear multiplied with higher powers of Za. 

The dependence of the (r 2 ) coefficient on the assumed proton charge distribution has been shown 
to be weak 0, \^\. This was demonstrated by integrating numerically the Dirac equation using 
various proton charge distributions: exponential, Gaussian and Yukawa (see Table HIj) . Another 
study integrating the Dirac equation numerically using the finite-size Coulomb potential given by 



the proton charge distribution extracted from the measured form factor (parametrized in |58y ) 



yields a total finite-size effect of [34 ] 



A-Efinite size = -5.19934r E 2 + 0.03458r E 3 + 0.00011r E 4 + 0.00027rl log (r E ) - 0.00005r4 log (r E ) . 

(18) 

Equation (I18D was attained by fitting the eigenvalues of the Dirac equation obtained for a finite- 
size Coulomb potential for various proton charge radii. Hence, it accounts for all-order finite-size 
effects. The choice of this parametrization was inspired by Eq. ([9]). 

The first coefficient of this equation is in agreement with b a + b c = —5.1973 — 0.00181 = 



—5.1991 meV/fm 2 of 28f] and the second one is compatible with Eq. (|15p . This implies that the 
two approaches, one starting from the Dirac equation with finite-size Coulomb potential, and the 
other one starting from the Schrodinger solution (with point-like Coulomb potential) complemented 
with relativistic and finite-size corrections are equivalent. The sum of the terms of Eq. f)18|) beyond 
r| and r| is only ~ 0.00004 meV, suggesting that the higher moments of the charge distribution 
(included in the form factor) do not affect significantly the prediction of the muonic hydrogen 
Lamb shift. 

When performing the integral of Eq. (|12p it resulted that the value is most sensitive to the 



cross sections measured at q 2 0.05 (GeV/c) 2 [5l|] and that there is also a considerable contribution 
from higher q 2 53J]. Thus an unexpected large value of R%-\ arising from contributions in the not 
measured low g 2 -region is tightly restricted by scattering data. The theoretical prediction of the 
finite-size terms given in Eq. Q may become more problematic if the higher moments of the 
charge distribution (r 4 ), (r 5 ), and (r 4 )(logr) etc. are large. Several of them have been evaluated 



using electron-scattering data |18j down to the lowest experimentally accessible exchanged photon 



momentum of g 2 ^ = 0.004 (GeV/c) 2 . The values reported in 18] for the higher moments are small 

? 3 



suggesting that the finite-size effect in muonic hydrogen is properly described by the (r 2 ) and Rf. 
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terms alone. 

However, extending arguments on the line of De Riijula, scattering experiments cannot totally 
exclude the existence of a "thorn" or a "lump" in the form factor Ge{q 2 ) at extremely low-g 2 



regime 59J] which could give rise to unexpectedly large higher moments of the charge distribution. 
Nevertheless the possibility to have such a behavior at low-q 2 is challenged by chiral perturbation 



theory (%PT) and vector meson dominance (VMD) models 24j] which account for the pion cloud 
in the low-energy regime. It was demonstrated by Pineda that is about 2-3 fm 3 in the leading 



chiral expansion term 27j, |54J . This sets tight constraints on the long tail of the proton charge 
distribution. Still it is important to further investigate the proton structure at very low q 2 with 
various techniques from xPT) to lattice QCD and VMD. 

Indelicate reevaluated the main radiative corrections including finite-size effects (see Table Hi]) . 
He computed the all-order finite-size corrections to the all-order one-loop eVP (one-loop eVP + eVP 
iteration) and to the Kallen-Sabry (including eVP iteration) term, using the charge distribution 



from the measured form factor parametrized by Arrington et al. [5Js|. Moreover Indelicate and 



Mohr calculated the finite-size correction to the muon self-energy perturbatively, and confirmed 



their result by using also an all-order finite-size approach 57j. Finite-size corrections to higher- 
order radiative contributions are negligible in our context. 

The higher-order finite-size correction given in Table ED for Borie is AE. B ™ = 



-0.000123 meV 



281 ] originating from the terms in Eq. © not proportional to (r 2 ) and i? 3 , 



(2)- 

ho, Indclicato 



For Indelicate the higher-order finite-size correction A^°! t ™ cato = 0.00001(10) meV (Eq. (115) 



m 



34j ) is the sum of the terms in Eq. (I18j) beyond r\ and r\ and the radiative finite-size corrections 
beyond the (r 2 ) terms for one- loop eVP, Kallen-Sabry and muon self-energy. 

From Table [H] we see that the finite size contributions are (OPE: one-photon exchange, TPE: 
two-photon exchange) 

A£ n h nitcsizc =A£ pe +A^ a p d E +A£ TPE +Ai^ P d E +A£ fi h n ° itcsizc meV 

= -5.1994r| -0.0275r| +0.0333(20) -0.00062 r| + 0.00018 +0.00001(10) meV 

(19) 

resulting in 

A^finite size = -5.2275(10) r\ + 0.0333(20) + 0.0002(1) meV (20) 

where the uncertainty of the r E coefficient accounts for the difference between the all-order and per- 
turbative approach and uncertainties/approximations related with the proton charge distribution. 
Including the contributions summarized in Table U we obtain the total Lamb shift prediction 

A££ h = 206.0336(15) - 5.2275(10) r 2 +0.0333(20) meV 

= 206.0669(25) - 5.2275(10) r| meV. 

This can be compared with the results in "Old-Nature" 

A£ L h = 206.0573(45) - 5.2262 r| +0.0347 r| meV 
= 206.0779(45) - 5.2262 r 2 meV 
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and in Borie-v6 

AE^ = 206.0592(60) - 5.2272 r| + 0.0365(18) r| meV 

= 206.0808(61) - 5.2272 r| meV . 

The difference in the constant terms between Eq. (|2ip and Eqs. (|22p and (j23|) is originating mainly 
from item #18 in Table Ull (0.013 meV) which was double counted in "Old-Nature". Another 
double counting in "Old-Nature" was related with #21 and #12. The uncertainty of the proton 
structure independent term in Eq. (1211) in given mainly by the uncertainties of the radiative-recoil 
correction #24 and uncalculated higher-order terms. 



IV. 2S-HFS AND THE ZEMACH RADIUS 

The interaction between the bound particle with the magnetic field induced by a nucleus mag- 
netic moment gives rise to shifts and splittings of the energy levels termed hyperfine effects. In 
classical electrodynamics the interaction between the magnetic moments fi p , and /x^ of proton and 



muon respectively is described by [35] 



^clas s sical = _2^.^ (y>) (24) 

where <5(r) is the delta-function in coordinate space. A similar Hamiltonian can be derived in 
quantum field theory from the one-photon exchange diagram. Using the Coulomb wave function 
thisgive rise in first order perturbation theory to an energy shift for muonic hydrogen nS-states 
of 



28] 



3 

3n 3 m /i m p v ' /v ' ^' 2 L v ' ' 2 



Erfs(F) = Ifj ™ r (! + «)(! + F(F + 1 



AE Feimi - 



(25) 

ir„,„ 3- 



F(F + 1)-- 



where A-Epa-mi — 22.8320 meVj28l] is the Fermi energy, m p the proton mass, F the total angular 
momentum, k and the proton and muon anomalous magnetic moments, respectively. The F = 1 
state is shifted by 1/4 x 22.8320 meV whereas the F = state by -3/4 x 22.8320 meV (see Fig. [I]). 
Equation ()25[) accounts for the sum of the terms (hi) and (h4) in Table [TTTl The Breit term (h2) 
corrects for relativistic and binding effects accounted in the Dirac-Coulomb wave function but 
excluded in the Schrodinger wave function. 

Table ITO1 also summarizes the corrections arising from QED, recoil, nuclear structure, hadronic 
and weak interaction effects. The structure dependent corrections, scaling as the reduced mass of 
the system, become large in pp compared with hydrogen. The largest correction is thus given by 
finite-size effects which, in the non-relativistic limit, is given by the well known Zemach term (h20) 

A^Zcmach = -AEpcrmi ■ 2{Za)m r r z (26) 

where rz is the Zemach radius defined as 

rz = / d 3 r / dV rp E (r)p M (r - r') (27) 
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with pm and pe being the normalized proton magnetization and charge distributions, respectively. 
The convolution between charge and magnetic distribution in rz is a consequence of the interac- 
tion of the proton spin distributed spatially (given by the magnetic form factor) with the spatial 
distribution of the muon spin which is described by the atomic muon wave function. The latter is 
slightly affected, notably at the origin, by the charge-finite-size effect and thus by pe- In a quan- 
tum field framework this contribution arises from two-photon exchange processes. Similar to the 
above discussed situation for the Lamb shift, in the intermediate virtual state the proton may be 
"on-shell" but also "off-shell". Hence polar izability contributions need to be accounted for (h22). 
This term has the largest uncertainty which arises from the uncertainty of the polarized structure 
functions g\ and g<i (measured in inelastic polarized electron-proton scattering) needed as an input 
to calculate this contribution. Contrarily to the Lamb shift, for the HFS the calculation of the 
two-photon exchange diagrams via Compton scattering and dispersion analysis do not require any 
subtraction term 



60] 



The leading recoil correction (h23) to the HFS is generated by the same two-photon exchange 
diagram and is of order (Za)(m/M)E-p erm i where £?Fermi is the Fermi energy without contribution 



of the muon anomalous magnetic moment 35]. 



Radiative corrections which are not accounted in the anomalous magnetic moment are listed 
separately in Tab. IIII1 The largest radiative correction is related with the distortion of the wave 
functions caused by the Uehling potential, accounted in (h5), (h7), (hll) and (h25). Other cor- 
rections account for modifications of the magnetic interaction caused by the eVP in one and two 
photon exchange (h8), (h9) and (hlO), and vertex corrections caused mainly by the muon self- 
energy (hl3) and (hl4). 



The main HFS contributions have been confirmed and refined by Indelicate [34| by numerical 
integration of the hyperfine Hamiltonian with Bohr-Weisskopf (magnetization distribution) cor- 
rection using Dirac wave functions. The latter have been calculated for Coulomb finite-size and 
Uehling potentials. All-order finite-size, relativistic and eVP effects are thus accounted in the wave 
function. This calculation is performed for various te and rz, assuming exponential charge and 
magnetization distributions. From these calculations the terms (h3), (h6), (h20), (h21) and (h25) 
are obtained, showing good agreement with the perturbative results. It is interesting to note that 
the HFS shows a small dependence on r E given by the term (h21). The small constant terms in 
(h21) and (h25) account for the sum of higher order terms of a polynomial expansion in r F and rz 



(r|, r^rz, r^rz, r-&r\, r|r| etc.) of the numerical results obtained in [341 ] . 
Summing up all contributions in Table IIIII we get 

AE$ FS = 22.9763(15) - 0.1621(10)r z + AE%£ S meV. (28) 



Inserting the value AE^ S = 0.0080(26) meV 65J results is 

A^ h FS = 22.9843(30) - 0.1621(10)r z meV. (29) 



The uncertainty for the first term in Eq. (|28p considers differences between results from various 
authors and uncalculated higher-order terms. 
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V. CONCLUSION 

We have presented an update of the theoretical prediction of the Lamb shift and the 2S-HFS 
splitting in muonic hydrogen. Double counting of a few higher-order terms in previous compilations 
has been eliminated. No large missing term or error have been found. Finite-size and one-loop 
eVP have been studied also by performing numerical integration of the Dirac equation including 



the finite-size Coulomb and the Uehling potential [34f | . This all-order calculation confirms the 
perturbative results, for both the Lamb shift and the HFS. Combining these predictions with the 
two measurements in muonic hydrogen [l| will result in the determination of the proton charge 
and Zemach radii. 
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TABLE I: All known radius-independent contributions to the Lamb shift in /ip from different authors, 
and the one we selected. Our selection considers the results from Indelicato indicated by Eq. (number) 
in • Values are in mcV. The entry # in the first column refers to Table 1 in "Old-Nature" supplementary 
informations Q. The "finite-size to relativistic recoil correction" (entry #18 in which depends on the 
proton structure has been shifted to Table [TTl together with the small terms #26 and #27, and the proton 
polarizability term #25. 



SE: self-energy, VP: vacuum polarization, LBL: light-by-light, Rel: relativistic, NR: non-relativistic. 



# 


Contribution 


Pachucki [3, 4J 


Old-Nature [2] 


Boric-v6 [28] 


This work 


Ref. 




1 


NR one-loop electron VP (eVP) 


205.0074 












2 


Relativistic (rel.) corr. (Brcit-Pauli) 


0.0169° 












3 


Rel. one-loop cVP 




205.0282 


205.0282 


205.02821 


Eq. (54) in [34] 




19 


I^gI recoil corr to gVP ^ ^ ol ^ ^ 


(inch in #2) 6 


— 0.0041 


— 0.0041 


—0 00?08 c 


[30 311 




4 


Non-rel. two-loop eVP (Kallcn-Sabry) 


1.5079 


1.5081 


1.5081 


1.50810 


Eq. (57) in [24] 




5 


One-loop eVP in 2-Coulomb lines o^(zjct)^ 


0.1509 


0.1509 


0.1507 


0.15102 0- 


Eq. (60) in [34J 




7 


eVP corr. to Kallcn-Sabry 


0.0023 


0.00223 


0.00223 


0.00215 e 


Eq. (62) in [34], [411 




6 


Non-rel. three-loop eVP 


0.0053 


0.00529 


0.00529 


0.00529 


[41,42J 




9 


Wichmann-Kroll, "1:3" LBL scatt. 




-0.00103 


-0.00102 


-0.00102 


Eq. (64) in [341. [431 




10 


Virtual Delbriick, "2:2" LBL scatt. 




0.00135 


0.00115 


0.00115 


[35, 43] 




new 


"3:1" LBL scattering 






-0.00102 


-0.00102 


[43] 




20 


Muon self-energy and /a VP 


-0.6677 


-0.66770 


-0.66788 


-0.66761 


Eq. (72) and (75) in 


[34] 


11 


Muon self-energy corr. to eVP a 2 (Za) 4 


-0.005(1) 


-0.00500 


-0.004924 1 


-0.00254 


Eq. (29a) in [39] 9 




12 


eVP loop in self-energy a 2 (Za) 4 


-0.001 


-0.00150 






[35. 44-461 




21 


Higher-order corr. to /j-SE and /^VP 




-0.00169 


-0.00171 h 


-0.00171 


Eq. (177) in [40J 




13 


Mixed electron and muon VP (eVP+/iVP) 




0.00007 


0.00007 


0.00007 


[35] 




new 


eVP and /^VP in two Coulomb lines 








0.00005 


Eq. (77) in [24] 




14 


Hadronic VP a(Za) l m r 


0.0113(3) 


0.01077(38) 


0.011(1) 


0.01121(44) 


[47-49J 




IF, 


Hadronic VP a(Za) 5 m r 




0.000047 




0.000047 


[48, 49J 




16 


Rad corr. to hadronic VP 




-0.000015 




-0.000015 


[48. 491 




17 


Recoil corr. 


0.0575 


0.05750 


0.0575 


0.05747 


Eqs. (81) and (82) in [34] 


22 


Rel. recoil corr. (Za) 5 


-0.045 


-0.04497 


-0.04497 


-0.04497 


Eq. (83) in [34J 




23 


Rel. recoil corr. (Za) G 


0.0003 


0.00030 




0.0002475 


Eq. (85) and Tab. II 


m [34] 


new 


Rad. (only eVP) recoil corr. a(Za) 5 








0.000136 


Eq. (64a) in [39] 




24 


Rad. recoil corr. a(Za) n (proton SE) 


-0.0099 


-0.00960 


-0.0100 


-0.01080(100) [251 ' 






Sum 


206.0312 


206.02915 


206.02862 


206.03339(109) 







"This value has been recalculated to be 0.018759 meV [3(J. 

'This correction is not necessary here because in #2 the Breit-Pauli contribution has been calculated using a 

Coulomb potential modified by eVP. 
"Difference between Eq. (6) and (4) in E^ ) (2P 1/2 - 2S 1/2 ) - £^(2P 1/2 - 2S 1/2 ) = 0.018759 - 0.020843 = 

= —0.002084 meV (see also Table IV). Using these corrected values, the various approaches are consistent. 

Pachucki becomes 205.0074 + 0.018759 = 205.0262 meV and Borie 205.0282 - 0.0020843 = 205.0261 meV. 
d All order calculation gives more complete results. 
C A11 order calculation gives more complete results. 
^In Appendix C incomplete. 

9 Eq. (27) in [39( includes contributions beyond the logarithmic term with modification of the Bethe logarithm 

to the Uehling potential. The factor 10/9 should be replaced by 5/6. 
h Borie includes wave function corrections calculated in [4l|]. The actual difference between "Old-Nature" 

and Borie-v6 is given by the inclusion of the Kallen-Sabry correction with muon loop. 
"This was calculated in the framework of NRQED. It is related with the definition of the proton radius. 
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TABLE II: Proton structure dependent contributions to the Lamb shift in /xp from different authors, and 
the one we selected. Values are in meV, (rf) in fm 2 . The entry # in the first column refers to Table 1 in 
"Old-Nature" supplementary informations [2|]. Entry # 18 is under debate. 
TPE: two photon exchange, VP: vacuum polarization, Rel: relativistic. 



# 


Contribution 


EB-12 [281 


SK-11 [311 


KP-99 [3, J 


PI-12 [341 


Carroll [371 


This work 




Non-rel. finite-size 


-5.1973{r 2 ) 


-5.1975(r 2 ) 


-5.1975(r 2 ) 










Rel. corr. to non-rel. finite size 


-0.0018{r 2 ) 




-0.0009 meV 










Rel. finite-size 
















form factor 








-5.1994(r 2 ) 




-5.1994(r 2 ) 




exponential 










-5.2001(r 2 ) 






Yukawa 










-5.2000(r 2 ) 






Gaussian 










-5.2001(r 2 ) 






Finite size corr. to one-loop eVP 


-0.0110{r 2 ) 


-0.0110(r 2 ) 


-0.010(r 2 ) 


-0.0281(r 2 ) 




-0.0281(r 2 ) 




Finite size corr. to one-loop eVP-iter. 


-0.0165{r 2 ) 


-0.0170(r 2 ) 


-0.017(r 2 ) 


(incl. in -0.0281) 








Finite-size corr. to Kallcn-Sabry 


b 






-0.0002(r 2 ) 




-0.0002(r 2 ) 


new 


Finite size corr. to muon self-energy 


(0.00699) " 






0.0008(r 2 ) 




0.0008(3)(r 2 ) << 




TPE finite-size contr. A_E T pe [231: 












0.0333(20) meV 




elastic (third Zemach) e 
















measured form factor 


0.0365(18)(r 2 ) 3/2 






0.0353(r 2 ) 3/2 f 




(incl. above) 




exponential 






0.0363(r 2 ) 3/2 




0.0353(r 2 ) 3/2 






Yukawa 










0.0378(r 2 ) 3/2 






Gaussian 










0.0323(r 2 ) 3/2 




25 


inelastic (polarizability) 


0.0129(5) meV [551 




0.012(2) meV 


0.0129(10) 




(incl. above) 


new 


Rad. corr. to third Zcmach, afZct) 5 


-0.00062(r 2 ) 










-0.00062(r 2 ) 


26 


eVP corr. to polarizability a(Zcx) 








0.00019 meV 




0.00019 meV [49J 


27 


Sclf-E corr. to polarizability a(Za) 5 








-0.00001 meV 




-0.00001 meV [491 


18 


Finite-size to rel. recoil corr. 


(0.013 meV) 3 




h 






(incl. in AEtpe) 




Higher-order finite-size corr. 


-0.000123 meV 






0.00001(10) moV 




0.00001(10) meV 




2P1/2 finite-size corr. 


-0.0000519(r 2 ) 1 






(incl. above) 


(incl. above) 


(incl. above) 


Corresponds to Eq. (6) in Q] which accounts only for the main terms in Frel 


and Fnrel- 







'This contribution has been accounted already in both the -0.0110 meV/fm 2 and -0.0165 meV/fm 2 coefficients. 
c Given only in Appendix C. 

d This uncertainty accounts for the difference between all-order in Za and perturbative approaches [5t1 |. 
'Corresponds to Eq. (Tl6)) . 

^This value is slightly different from Eq. (|18[) because here an all-order in finite-size AND an all-order in eVP 

approach was usecl^ 
9 See Appendix F of [50]. This term is under debate. 

''Included in A£ T pe. This correction of 0.018 - 0.021 = -0.003 meV is given by Eq. (64) in 0] and Eq. (25) in 0]. 

This correction is also discussed in [i^ where the 6/7 factor results from 0.018/0.021. 
l Eq. (6a) in 0. 
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TABLE III: All known contributions to the 2S-HFS in /zp from different authors, and the one we selected. 
Values are in meV, radii in fm. Our selection considers the results from P. Indelicato ■ 





Contribution 


Martynenko[32l 


Boric- v6 [28] 


This work 


Rof. 




hi 
h2 
h3 


Fermi energy, (Za) 4 
Breit correction, (Za) G 

Dirac energy (includes Breit corr. in all-order) 


22.8054 
0.0026 


22.8054 
0.00258 


22.807995 


Eq. (106) 


in [34] 


h4 


[A anom. magnetic moment corr., a(Za) 4 , a(Za) 4 


0.0266 


0.02659 


0.02659 






h5 
h(i 
h7 


eVP in second-order pert, theory, a(Za) 5 (eyp2) 
All-order eVP corr. 

Two-loop correction to Fermi-energy (eyp2) 


0.0746 


0.07443 
0.00056 


0.07437 
0.00056 


Eq. (108) 


in [34] 


h8 
h!) 
hlO 


One-loop eVP in I7 interaction, a(Za) 4 (eypi) 
Two-loop eVP in I7 interaction, a 2 (Za) 4 (evpi) 
Further two-loop eVP corrections 


0.0482 
0.0003 


0.04818 
0.00037 
0.00037 


0.04818 
0.00037 
0.00037 


[61, 62] 




hll 
hl2 


fiVP (similar to eyp2) 
fiVP (similar to evpi) 


0.0004 


0.00091 
(inch in hl3) 


0.00091 
(incl. in hl3) 






hl3 
hl4 


V CI I CX, OiyZj a j 

Higher order corr. of (hl3), only part with ln(a) 




-U.UUOll 

-0.00017 


-U.UUOll 

-0.00017 


[63] 




hl5 
hl6 
hl7 


fi self-energy corr. with proton structure, a(Za) 5 
Vertex corrections with proton structure, a(Za) 5 
"Jellyfish" diagram corr. with p structure, a(Za)^ 


0.0010 
-0.0018 
0.0005 










hl8 
hl9 


Wndrnn VP rv Q 
nduiuit v Jr , ex 

Weak interaction contribution 


n nnn^c 1 ^ 

U.UUUO^ 1 ) 

0.0003 


n nnnfifV i n\ 

U . UUUDU^ 1UI 

0.00027 


n nnnfin/' i n\ 

U . UUUDU^ 1UI 

0.00027 


[64] 




h20 
h21 


Finite-size (Zemach) corr. to ASpcrmi, (^ a ) 5 
Higher-order finite-size corr. to A£^p erm i 


-0.1518 6 


-0. 16037 r z 


-0. 16034 r z 
-0.0022 r E 2 +0. 0009 


Eq. (106) 
Eq. (106) 


in [34] 
in [34] 


h22 


Proton polarizability, (Za) 


0.0105(18) 


0.0080(26) 


0.00801(260) 


[65. 661 




h23 


Recoil corrections 




U. Ul J. ZiO 


0.02123 


[60] 




h24 
h25 
h26 


eVP contribution -f proton structure corr., a G 
eVP corr. to finite-size (similar to eyp2) 
eVP corr. to finite-size (similar to eypi) 


-0.0026 


-0.00114 
-0.00114 


-0.0018 r z -0. 0001 
-0.00114(20) 


Eq. (108) 


in [34J 


h27 
h28 


Proton structure corr., a(Za) J 

Rel.+Rad recoil corr. with proton AMM, a 6 


-0.0017 c 
0.0018 












Sum 


22.8148(20) d 


22.9839(26) 
0.1604 r z 


22.9858(26) 

-0.1621(10) r z -0.0022(5) r| 






Sum with r E = 0.841 fm and r z = 1.045 fm [181 


22.8148 meV 


22.8163 meV 


22.8149 meV 







"Includes a correction a(Za) 5 due to //VP. 
'Calculated using the Simon et al. form factor. 
c Is this term included in (h26)? 

d The uncertainty is 0.0078 meV if the uncertainty of the Zemach term (h20) is included (see Table II of [j^h 



